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Abstract 

The wave function at the origin (WFO) is an important quantity 
in studying many physical problems concerning heavy quarkonia. 
However, when one used the variational method with fewer pa¬ 
rameters, in general, the deviation of resultant WFO from the 
’’accurate” solution was not well estimated. In this paper, we dis¬ 
cuss this issue by employing several potential forms and trial wave 
functions in detail and study the relation between WFO and the 
reduced mass. 
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I. Introduction 


Recently, the wave function at the origin for the S-wave bound state of a heavy quark- 
antiquark system once again attracts physicists’ attentions0’ Hi. This is because that it is not 
only a very important quantity for calculating spin state hyperfine splitting, but also crucial 
to evaluating the production and decay amplitude of the heavy quarkonium. Within the 
context of the non-relativistic potential model, Refs. and demonstrated the numerical 
results of WFO of the S-wave cc , be and bh systems and compared with those obtained in 
various ’’successful” potential models. 

As well known, except the Coulomb and the harmonic oscillator potentials, there are 
few potentials which bound state problems can be analytically solved. For solving these 
non-analytically soluble bound state problems, one has to use approximations. Numerically 
solving Schrodinger equation is the most powerful method which can reach most required 
accuracy. But the numerical method has some defects, for instance, it cannot give analytical 
expressions for further discussion. Moreover, all numerical method for the central potential 
is only available for the v{r) which has the singularity less than ^ when r approximate to 
0, therefore, it definitely fails as ^ exist in the potential. R is unfortunately, the case is 
encountered usually in calculating the fine-splitting of the P state. 

The perturbative method is another approximation method which has most extensively 
been used. However, the practical applicability of the perturbative expansion in many cases 
is limited due to divergence, and the ranges of perturbation parameter values are usually 
restricted by the convergence requirements. Moreover, in the perturbation method, the 
treatment for wave functions is much more difficult than that for energy eigenvalues. As 
it was declared by MacClary and Byers that it is not simple to obtain the wave function 
correct to the order of because the perturbative correction to the wave function should 

be given by an infinite sum over all states of the nonrelativistic Hamiltonian. As an example, 
in the positronium case, this sum diverges because of the divergence of the potential at the 
origin. Besides, for the singular ^ potential, the perturbative method is the one used most 
frequently yet. The splitting of the energy level is obtained by calculating an integral. How¬ 
ever, as mentioned in Ref.(5|, such a singular potential would lead to an exotic result which 
does not correspond to the real physics. Therefore, the meaning of the result given by the 
perturbative calculation is not evident actually. 

Of course, as a physical state which can be measured in the experiment, there should 
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be a theoretically derived solution to match, as long as the employed model is correct. The 
problem is how to get it. For instance, the energy levels of the triplet P states of the cc 
system must have definite measured values. In the theoretical calculation, due to the non¬ 
soluble nonperturbative QCD effect, one has to rely on specific models. The non-relativistic 
quark potential model is one of them. Under the tree diagram approximation and the non- 
relativistic approximation, a potential term of r~^ appears. This is simply because that even 
in the perturbative framework, one does not collect (actuarially is not able to do so) all 
the high order diagrams in the perturbative expansion and all the high order terms in the 
non-relativistic reduction. Then the question is how to catch the major character of the real 
physics. 

There were lots of conscientious attempts have been made. It is noteworthy that in 
Gupta’s papers® 0, the non-relativistic reduction was performed with respect to the power 
of /E‘^ so that the singular r~^ potential term can be avoided. Our recent study indicates 
that the variational method can give some interesting sight to the singular ^ problem. 

As well known that the variational methods a widely employed approach. In principle, 
by using this method one can leave the above mentioned problem alone and get the major 
content of the real physics. We will discuss this issue in our successive paper Q. In spite of 
this, the variational method has more advantages. It can give an analytical expression of the 
wave function. In particular, if there is only a single parameter in the trial wave function, 
the resultant wave function has a simple form. Then it is very convenient in the practical 
application and physical discussion. 

The variational method has extensively been used and seems to be successful in many 
aspects. Although most of works paid their attentions on seeking out accurate energy eigen¬ 
values and seldom discussed the wave function, in particular WFO, Ref. Q simply discussed 
the wave function by employing the multi-Gaussian trial wave function. 

In many simple cases, the Gaussian-type or the exponential-type functions with a single¬ 
parameter were taken as the trial wave function to discuss the ground or excited states. In 
general, the accuracy of the resultant energy eigenvalues is satisfactory. However, by looking 
at the wave function, one would find that although the resultant wave function sometimes 
may not deviate much from the ’’accurate” solution (usually it can be obtained by solving 
the Schrodinger equation numerically) in the long-range part, but this deviation turns larger 
and larger when r ^ 0. 
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In this work, we would study this problem in detail with various potentials which were 
used frequently in studying heavy quarkonia and put more emphasis on estimating WFO. 
The adopted trial wave functions were used in literatures except the last one. 

The paper is organized as following. After the introduction, the variational method by 
using the single—parameter trial wave function is discussed in Sec.II. In Sec.Ill, variational 
calculation by employing multi-parameter trial wave functions is presented. The variational 
study on the 2S state of the cc system is further shown in Sec.IV, and in Sec.V, the relation 
between WFO and the reduced mass for various trial wave functions is investigated. In the 
last section, the discussion and conclusion are given. 

II. Variational method with single—parameter trial wave function 


There are many potential models which can fit the experimental spectra of the heavy 
quarkonia with certain accuracy. In the rest of the paper, we call them as ’’successful” 
potential models. 


Within the framework of the non-relativistic potential model, the S-state wave function 
'0(r) of the heavy quarkonium satisfies the Schrodinger equation 

H-ipir) =——Aijjir)+ V{r)il^(r) = (1) 

2 // 

where H is the Hamiltonian of the quarkonium, V{r) denotes the central potential between 
quark and antiquark, E represents the energy eigenvalue, and /r is the reduced mass. In gen¬ 
eral, one can numerically solve the Schrodinger equation to obtain E and 'ipi'r’) simultaneously, 
and calculates ■0(0) (WFO) in terms of the average value of ^ ||3|, i.e. 


9 u dV 
= — < — > 
27r dr 


( 2 ) 


The values of the squared WFO for various ’’successful” potential models were listed in Refs, 
in, 1^. (note: |R(0)p is equal to 47r|'!/;(0)p in the Table 1 of Ref.ll|. )• 


To solve Eq.(l) by using the variational method, one needs to choose a suitable trial wave 
function V’(^; c) with N independent parameters {c} = {ci, C 2 , ■ ■ ■, cw} first and then to seek 
out a set of parameters {cq} = {cio,i = 1,2, ■ • ■,N} which minimizes the expectation value 
of Hamiltonian, namely 


E{c) =< H >= 


< V’(c)|H~|V:(c) > 
< V’(c)|V’(c) > 


(3) 
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The minimum value E{cq) gives an upper limit of the ground state energy. In the use of the 
variational method, one wishes to obtain E{cq) which is as close as possible to the ’’accurate” 
solution with the minimum number of parameters. 

In this work, what we concern is how close to the ’’accurate” solution the resultant wave 
function can be, when E{co) is satisfactorily close to the ’’accurate” energy eigenvalue. In 
particular, we would try to find the accuracy of WFO for various trial wave functions, which 
obviously affects the application of the variational method. 

In this paper, we choose three most popular and ’’successful” models listed in Refs.|j|, ||], 
so that the conclusion could be more general. These models are |^: 

(1) . Cornell potential IT^]: 

V{r) = -\— + kr, (4) 

6 r 

with as = 0.39, k = 1/2.34^(GeF)^ and the mass of c quark rUc = l.MGeV. 

(2) . Martin potential |pl|: 

V{r) = kr^'^, (5) 


with k = 6.898 and rric = l.SGeV. 

(3). Logarithmic potential [p!^: 

V (r) = klog{r), 


( 6 ) 


with k = 0.733 and rric = 1.5GeV. 

In this section, we choose the simplest trial wave function in which there is only one 
variational parameter to study the 15 state of cc. The general form of such trial wave 
function is written as 


Ariai{r) = N e “ 


(7) 


where 


N = 


b (2a )b 

4vrr(f) 


1 

2 . 


( 8 ) 


■^Indeed, k in Eqs.(5) and (6) should have proper dimensions, but in our work, it is not important because 
they would be attributed into the normalization factors. 
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is the normalization constant, a denotes the variational parameter which will be fixed by min¬ 
imizing the expectation value of Hamiltonian and b is the model parameter which determines 
the type of the trial wave function. In practice, we select following four trial wave functions: 

(1) . 6 = 1, namely N ^ (hydrogen wave function or exponential wave function). It is 
the solution of the Coulomb potential model. 

(2) . 6 = 2, namely N ^ (harmonic oscillator wave function or a Gaussian wave 
function). 

(3) . 6 = |, namely N e~°‘ This function was used by Gupta 

(4) . 6 = |, namely N e~°‘ , This is a newly proposed trial wave function, and we will 

pay more attention on it. 


In order to fully understand the accuracies of the variational results in three different 
potential cases, we calculate four quantities by using four different trial wave functions and 
the corresponding ’’accurate” results by solving the Schrodinger equation numerically as well. 
These quantities are: 

(1) . the energy eigenvalue E (note: we do not try to fit the experimental spectrum here, 
because it is not the aim of this work). 

(2) . average radius < r >. 

(3) . average value of the inversed radius < ^ >. 

(4) the squared WFO |V'(0)p. 


For each calculated quantity q we give a relative deviation 6q, which is defined as 

(,Qvar Qtrue) 


6q = 


Qtrue 


where q^uar is the variational result and qtrue the ’’accurate” value. 


In the case of the Gornell potential (4), kinetic energy < T > and potential energy <V> 

are: 


< T > 

< F > 


(2a)f62F(2 + i) 

8/^r(f) ’ 

-4(2a)f«,r(|)+3fcF(|) 

3(2a)ir(f) 


(9) 

( 10 ) 


Then we can obtain the expectation value of Hamiltonian and consequently an algebraic 
equation, which is used to determine a, 

[36^r(2 ^)]x^ - [16asfjT{^)]x‘^ - 12/c^F(^) = 0, 
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where x = (2a) t. It is very easy to solve this equation. If we rewrite it in the following form: 

+ A.2X^ + Aq = 0 , ( 11 ) 

the real solution of x can be expressed as: 


A 2 23 An 

^real TTU 1“ TTZ “1“ 


B 


SAs 3A3B 3 25743’ 

where B = {Bq + —AA^ + Bq)3 and Bq = —2^2 “ 277 I 0 A 3 . 

In the case of Martin potential (5), the potential energy reads: 

^r(¥) 

< F >= ^ ^ ^ 


(2a)-r(|) 

Therefore, the equation for determining a is quite simple. The solution is 


a = 


1 0.4^fcr(M) 

2^ 62^(2+1) J 


( 12 ) 


(13) 


(14) 


In the case of logarithmic potential (6), by the similar procedure we obtain 

^^h±jL[AL,i 

ry r (2 + i)' ^ 

All the numerical results are listed in Tables 1.1 to 1.3. In order to make comparison, we 
write the corresponding ’’accurate” results in the table captions and list relative deviations 
in the corresponding tables. 



Table 1.1 

The variational results with a single parameter trial function in the Cornell potential case. The values 
listed in parentheses are relative deviations. The ’’accurate” results are: E{1S) = 0.257526 GeV, < r > = 
1.7073 GeV-^, < 7 >= 0.80848 GeV and = 0.116054 GeV^. 


b 

E 

< r > 

< 

7 > 

iV’(o)r 

1 

0.272795 ( 0.059) 

1.7908 ( 0.049) 

0.8376 

( 0.036) 

0.197052 ( 0.61) 

2 

0.280039 ( 0.087) 

1.7234 ( 0.009) 

0.7388 

( -0.098) 

0.050403 ( -0.57) 

'6 

9 . 

0.259785 ( 0.0088) 

1.6989 (0.0049) 

1 0.7908 ( -0.022) 

0.082911( -0.29) 

4 

3 

0.257809 (0.0011) 

1.7083 ( 0.0006) 

1 0.8083 

(- 0 . 0002 ) 

0.103334( -0.11) 


Table 1.2 


The variational results with a single parameter trial function in the Martin potential case. The values 
listed in parentheses are relative deviations. The ’’accurate” solutions are: E{1S) = 7.5605 GeV, < r > = 
1.72332 GeF“\ < 7 >= 0.782946 GeV, \'ip{0)\^ = 0.0778779 GeV^. 
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b 

E 

< r > 

< 7 > 

\mv 

1 

7.5871 ( 0.0035) 

1.8606 (0.079) 

0.8064 ( 0.030) 

0.166933 ( 1.14) 

2 

7.57509 ( 0.0019) 

1.7151 (-0.0048) 

0.7423 (-0.052) 

0.051445 (-0.34) 

d 

7.5607 (3>k 10"U 

1.7194 (-0.0023) 

0.7814 (-0.0020) 

0.079984(-0.027) 

3 

7.5622 (0.0002) 

1.7415 (0.011) 

0.7928 (0.013) 

0.097596 (-0.25) 


Table 1.3 


The variational results with a single parameter trial function in the logarithmic potential case. The 
values listed in parentheses are relative deviations. The ’’accurate” solutions are: E{1S) = 0.730733 GeV, 
<r>= 1.87535 GeV~\ < 7 >= 0.723538 GeV, \'tp{0)\^ = 0.0633063 GeV^. 


b 

E 

< r > 

< 7 > 

iv>(o)r 

1 

0.754098 (0.032) 

2.0231 ( 0.079) 

0.7414 (0.025) 

0.129748 (1.05) 

2 

0.747750 (0.024) 

1.8639 (-0.0061) 

0.6831 (-0.056) 

0.039846 (-0.37) 

d 

‘2 

0.731223 (0.0007) 

1.8668 (0.0046) 

0.7197 (-0.0053) 

0.062492(-0.013) 

4 

3 

0.731726 (0.0014) 

1.8910 (0.0084) 

0.7301 (0.0091) 

0.076171(-0.20) 


III. Trial functions with two, three and four parameters 

The trial wave function with two or more variational parameters may have various forms. 
The most straightforward one is to multiply Eq.(7) by a polynomial of r, namely 

•0(r) = (co + Cl r + C 2 + ••• + c^ r”) e““(16) 

where a and ci, C 2 , • • •, are variational parameters, and cq can be fixed by the normalization 
condition. Using the standard procedure discussed above, we calculate the above mentioned 
four quantities in term of the trial wave functions with two, three and four parameters in 
the three potential cases and four b value cases. The numerical results are listed in Tables 
2.1—2.3, 3.1—3.3 and 4.1—4.3, respectively. 


Table 2.1 


The variational results with the two-parameter trial wave function in the Cornell potential case. The 
values listed in parentheses are relative deviations. The ’’accurate” results are: E{1S) = 0.257526 GeV, 
<r>= 1.7073 GeV~\ < 7 >= 0.80848 GeV and I'iAlO)!^ = 0.116054 GeV^ 


b 

E 

<r > 

< 7 > 

\mv 

1 

0.265416 (0.031) 

1.7562 (0.029) 

0.8296 (0.026) 

0.158048 ( 0.36) 

2 

0.271382 (0.054) 

1.7098 (0.0015) 

0.7585 (-0.062) 

0.065156 (-0.44) 

3 

2 

0.258738 (0.0047) 

1.6991 (0.0043) 

0.7959 (-0.016) 

0.091133(-0.22) 

4 

3 

0.257624 (0.0004) 

1.7080 (0.0004) 

0.8088 (0.0003) 

0.108602(-0.064) 


Table 2.2 
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The variational results with the two-parameter trial wave function in the Martin potential case. The 
values listed in parentheses are relative deviations. The ’’accurate” solutions are: E{1S) = 7.5605 GeV, 
<r>= 1.72332 GeF■^ < 7 >= 0.782946 GeV, \'ip{0)\^ = 0.0778779 GeV^. 


b 

E 

< r > 

< 7 > 

i^(o)r 

1 

7.56173 (0.0001) 

1.7337 (0.0060) 

0.7847 (0.0022) 

0.069440 (-0.11) 

2 

7.56873 (0.0011) 

1.7110 (-0.0007) 

0.7569 (-0.033) 

0.064430 (-0.17) 

l-i 

9. 

7.56068 (2 * 10“”) 

1.7213 (-0.0012) 

0.7830 (4 lO""*) 

0.082579 (0.060) 

4 

3 

7.56054 (5 * 10"*’) 

1.7222 (-0.0006) 

0.7832 (0.0001) 

0.0785065 (-0.0081) 


Table 2.3 


The variational results with the two-parameter trial wave function in the logarithmic potential case. The 
values listed in parentheses are relative deviations. The ’’accurate” solutions are: E{1S) = 0.730733 GeV, 
<r>= 1.87535 GeT■^ < 7 >= 0.723538 GeV, \'ip{0)\^ = 0.0633063 GeT®. 


b 

E 

< r > 

< 7 > 

mor 

1 

0.731264 (0.0007) 

1.88523 (0.0053) 

0.7246 (0.0015) 

0.058204 (-0.081) 

2 

0.740704 (0.014) 

1.8587 (-0.0089) 

0.6968 (-0.037) 

0.050334 (-0.21) 

3 

2 

0.731022 (0.0004) 

1.8707 b0.0025) 

0.7233 (-0.0017) 

0.066473(0.050) 

4 

3 

0.730798 (9 * 10-”) 

1.8741 b0.0007) 

0.7236 (0.0001) 

0.064889(0.025) 


Table 3.1 


The variational results with the three-parameter trial wave function in the Cornell potential case. The 
values listed in parentheses are relative deviations. The ’’accurate” results are: E{1S) = 0.257526 GeV, 
<r>= 1.7073 GeV-^, < 7 >= 0.80848 GeT and |iA(0)P = 0.116054 GeT^ 


b 

E 

<r > 

< 7 > 

iv>(o)r 

1 

0.257834 (0.0012) 

1.7096 (0.0013) 

0.8096 (0.0014) 

0.108931 (-0.061) 

2 

0.257908 (0.0015) 

1.7063 (-0.0006) 

0.8059 (-0.0032) 

0.105760 (-0.089) 

3 

2 

0.257665 (0.0005) 

1.7073 (-4>k lO-") 

0.8078 (-0.0008) 

0.105696 (-0.089) 

4 

3 

0.257623 (0.0004) 

1.7079 (0.0004) 

0.8086 ( 0 . 0001 ) 

0.108278(-0.067) 


Table 3.2 


The variational results with the three-parameter trial wave function in the Martin potential case. The 
values listed in parentheses are relative deviations. The ’’accurate” solutions are: E{1S) = 7.5605 GeV, 
<r >= 1.72332 GeT’b < 7 >= 0.782946 GeT, |^/>(0)p = 0.0778779 GeT®. 


b 

E 

< r > 

< 7 > 

mor 

1 

7.56082 (4* 10"”) 

1.7264 (0.0008) 

0.7843 (0.0017) 

0.070081 (-0.10) 

2 

7.56077 (4* 10~”) 

1.7230 (-0.0002) 

0.7848 (-0.0023) 

0.090033 (0.16) 

3 

2 

7.56050 ( 6 * 10~') 

1.7233 (5>k lO-®) 

0.7829 (-3>k 10~”) 

0.077090 (-0.010) 

4 

3 

7.56069 (2* 10~”) 

1.7236 (-0.0002) 

0.7818 (0.0014) 

0.081307(-0.044) 







Table 3.3 


The variational results with the three-parameter trial wave function in the logarithmic potential case. The 
values listed in parentheses are relative deviations. The ’’accurate” solutions are: E{1S) = 0.730733 GeV, 
<r>= 1.87535 GeV~^, < 7 >= 0.723538 GeV, \'tp{0)\^ = 0.0633063 GeV^. 


b 

E 

< r > 

<7> 


1 

0.730759 (4* 10"'’) 

1.8762 (0.0004) 

0.7236 ( 8 >k lO-'’) 

0.065321 (0.032) 

a 

0.730964 (0.0003) 

1.8757 (0.0002) 

0.7245 (-0.0014) 

0.072312 (0.14) 

H 

0.73088 (0.0002) 

1.8745 (-0.0004) 

0.7240 (0.0006) 

0.0676953(0.069) 


0.730796 (9* 10~'’) 

1.8743 (-0.0006) 

0.7238 (0.0004) 

0.0651086 


Table 4.1 


The variational results with the four-parameter trial function in the Cornell potential case. The values 
listed in parentheses are relative deviations. The ’’accurate” results are: E{1S) = 0.257526 GeV, < r > = 
1.7073 GeV-^, < 7 >= 0.80848 GeV and |iA(0)P = 0.116054 GeV^. 


b 

E 

< r > 

< 7 > 

i^(o)r 

1 

0.257599 (0.0003) 

1.7093 (0.0012) 

0.8086 ( 0 . 0002 ) 

0.114842 (-0.010) 

2 

0.257667 (0.0005) 

1.7021 (0.0003) 

0.8096 (-0.0014) 

0.111811 b0.037) 

‘2 

0.257566 (0.0005) 

1.7068 (-0.0003) 

0.8077 bO.OOlO) 

0.106024(-0.086) 

4 

3 

0.257581 (0.0002) 

1.7074 ( 6 * 10”'’) 

0.8085 (5* lO””) 

0.110470 (-0.048) 


Table 4.2 


The variational results with the four-parameter trial wave function in the Martin potential case. The 
values listed in parentheses are relative deviations. The ’’accurate” solutions are: E{1S) = 7.5605 GeV, 
<r>= 1.72332 GeT■^ < 7 >= 0.782946 GeV, \'ip{0)\^ = 0.0778779 GeV^. 


6 

E 

< r > 

< 7 > 

\mv 

1 

7.56051 (1 10”'’) 

1.7233 (-2 * 10””) 

0.7830 (2 * 10””) 

0.0791203(0.016) 


7.56068 (2>k 10””) 

1.7192 (-0.0024) 

0.7852 (0.0029) 

0.0860623 (0.11) 

m 

7.56051 (1 10””) 

1.7226 (-0.0004) 

0.7830 ( 8 * 10””) 

0.07723 (0.0083) 

B 

7.56051 (1 10”*’) 

1.7232 (- 8 * 10””) 

0.7831 (0.0001) 

0.0779022 (0.0003) 


Table 4.3 


The variational results with the four-parameter trial wave function in the logarithmic potential case. The 
values listed in parentheses are relative deviations. The ’’accurate” solutions are: E{1S) = 0.730733 GeV, 
<r>= 1.87535 GeT■^ < 7 >= 0.723538 GeV, \'ip{0)\^ = 0.0633063 GeV^. 


B 

E 

< r > 

< 7 > 

mor 

1 

0.730737 (6 * 10””) 

1.87626 (0.0005) 

0.7234 (-0.0001) 

0.063531 (-0.0036) 

B 

0.730913 (0.0002) 

1.87499 (-0.0002) 

0.72366 (-0.0002) 

0.069626 (-0.10) 

B 

0.730738 (7* 10””) 

1.8754 (4* 10””) 

0.7236 (0.0001) 

0.064408 (0.017) 

B 

0.730732 (1 * 10””) 

1.8754 (2* 10””) 

0.7236 (0.0004) 

0.0629413 (-0.0058) 
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Another possible type of trial wave function with multi-parameters can be chosen in the 
form of the superposition of two or more above mentioned trial wave functions with different 
6 or a values, respectively. We give an example in the following. From the resultant |'0(O)P 
in Tables 1.1—1.3, it is easy to see that in the exponential-type trial wave function (b = 1) 
case, |^/;(0)p values are larger than ’’accurate” ones (the deviation is positive), while in the 
other three trial wave function cases (b = 2,3/2 and 4/3), the corresponding results are less 
than ’’accurate” ones (negative deviation). Therefore, one can compromise the deviations by 
combining two different types of trial wave functions with which the positive and negative 
deviations appear, respectively. After testing various combinations, we find that the results 
are similar to those by using three parameters in the last section. In fact, this new trial 
wave function has three variational parameters too. In the following, we list two specific 
combinations and corresponding results. 

The first example is the mixture of a exponential-type wave function and a Gaussian wave 
function, i.e. 

2 2 

R{r) = Cae~^ -h Coe~^^^. (17) 

For the Cornell potential, the relative deviation of energy is 0.002 and the relative devia¬ 
tion of squared WFO is 0.099. 

The second example is the mixture of an exponential-type function and a Gupta’s function, 

i.e. 

i?(r) = Cae“« -h Cg . (18) 

For the Cornell potential the relative deviation of energy is 0.00015 and the relative 
deviation of squared WFO is 0.025. 

Apparently, with the same number of variational parameters, these combined trial wave 
function can give better description on both energy and WFO. 

IV. 25 state of cc 

Based on the results of the 15 state, it is easy to discuss the 25 state. We first select 
a normalized trial wave function which is orthogonal to the 15 wave function. It can be a 
single-parameter or multi-parameter function. But we find that it is quite difficult to obtain 
a highly accurate WFO. In order to make sense, we demonstrate two examples. 
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The first example is that we take the 15 trial wave function as 

Ris{r) = N (19) 

where N is the normalization constant and a is determined by the variational method. Then 
we choose the 25 trial wave function to be 

R 2 s{r) = (co + Cl r + C 2 r^) , (20) 

where a takes the same value as that in Ris{r). By considering the orthonormal condition, 
only one parameter remains free. This parameter can be fixed by the variational method. 
The obtained results show that the relative deviation of energy is 0.004, while the relative 
deviation of the squared WFO of the 2S state is 0.25. 

Moreover, if we take 

R 2 s{r) = (co + cir + C 2 r^ + csr^) , (21) 

there are two variational parameters in the trial wave function. The resultant relative devi¬ 
ations of energy and squared WFO turn to be 0.0009 and 0.11, respectively. 

Similar to the 15 state trial wave function, if we take a trial wave function with four 
parameters, the resultant relative deviation of the squared WFO of the 25 state is less than 
0.05. 


V. Relationship between WFO and reduced mass 


There were lots of discussions concerning the relation between WFO and the reduced mass 
in the heavy quarkonium system [1,2]. Because the variational method can give the analytical 
expression of the wave function, consequently the exact value of WFO, the approximate 
relation between WFO and the reduced mass can be obtained. For example, if the trial 
wave function with parameter a for the 15 state, Eq.(7), is chosen, the squared WFO can be 
written as 

In the Martin potential case, substituting (14) into (22), we obtain 


.0.4A:F( 


M' 

b - 


47rF(|)^62r(2+i) 


/i]^. 


(23) 
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This relation is similar to that given by the simple scaling arguments for the power-law 
potential, namely 


2 



(24) 


When v = 0.1, Eqs.(23) and (24) coincide with each other. It clearly shows that although the 
variational method is an approximation, it retains the main characteristics of the solution. 
This observation encourages us to apply the variational method to the estimation of WFO. 


As shown in Table 1.1, when 6 = |, the squared WFO of the IS state of cc has the least 
deviation(0.027). By taking rri}, = 5.174GeF [1T|, one obtains 0.361475 for the squared WFO 
of the IS" state of bb. In the Be meson case, the calculated reduced mass is 1.3336 GeV, and 
the squared WFO of the 15 state is 0.140547. The resultant relative deviations of WFO for 
both bb and are 0.027 which is the same as that for cc. 


In the logarithmic potential case, the similar treatment leads to 


4vrr(^ 


r 4fcr(f) 


/r]i. 


(25) 


It is easy to see that the relation between WFO and the reduced mass /r is consistent with 
that obtained from the scaling arguments for the power-law potential with i/ = 0. 


Again in the single-parameter trial wave function with ^ = | case, one finds the least 
deviation. By taking rric = 1.5GeV and mi, = 4.906GeF |^], the squared WFOs of the 
15 states of J/V', Tl' and Be are 0.062492, 0.36964 and 0.118462, and the corresponding 
relative deviations are about 0.013, respectively. This is a quite satisfactory and encouraging 
observation. 


In the Cornell potential case, by employing the single-parameter trial wave function with 
a given in Eq.(12), one can carry out a similar discussion. The obtained relation between the 
squared WFO and the reduced mass is much more complicated than that in Fqs.(23) and 
(25). We omit the detailed expression in the text but give following two interesting results, 
which are computed by taking ^ = | and the relative deviation within 10%. 

(1) . The curve of |V'(0)|^ as the function of the reduced mass ^ has the similar behavior 
as those in the former two cases. When jj, = 0.2 — 3.0GeV , the function /(/i) ~ can fit 
the curve very well. This corresponds to = —0.65 in Fq.(24). 

(2) . If we take me = 1.84 GeV and mh = 5.17 GeV [^, the conjecture of Fq. (8) in 
Ref.i 
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|V’fe5(0)P ~ |V'c5(0)|^-^ IV'fe5(0)|°-^ 

holds within 2.5%. 

VI. Conclusion and discussion 

In this paper, we carefully studied the variational method, especially in determining the 
binding energies, wave functions at the origin, average radii and etc. of the quarkonium. 
Retaining generality as much as possible, we employ several ’’successful” potential models 
to analyze. By comparing the numerical results obtained in terms of the variational method 
with those by solving Schrodinger equation numerically, we present the relative deviations by 
employing different trial wave functions with single- or multi-parameters in various potential 
models. 

In the variational method with fewer parameters, even the calculated binding energies and 
some transition matrix elements are accurate enough, it is not easy to obtain a very accurate 
WFO. Because of the phenomenological requirements, a certain accuracy is requested. Thus, 
one of our aims is to seek for a possible and simpler way to solve this problem. 

We first study the single-parameter trial wave function case. The results shown in Tables 
1.1-1.3 indicate that for the Cornell potential, the trial wave function with & = |, i-e. = 

4 

N e~°‘ can give the least relative deviation of squared WFO. The value of the deviation 
is about 0.11 which is not accurate enough, although the relative deviation of energy already 
reaches 10“^. For the Martin and the logarithmic potentials, the situations are better. When 
6 = |, one obtains the least values of 0.027 and 0.013 for the relative deviations of squared 
WFO, respectively, while the corresponding relative deviations of energy reach 10“^ and 
10“®, respectively. 

The accuracy of variational results can be improved when the number of the variational 
parameters are increased. This is true in general. If one adopts more than four variational 
parameters, he may expect very high accuracy for energy, but not always for WFO. It is much 
more difficult to further improve the accuracy of WFO than that of energy. The resultant 
accuracy of WFO seriously depends on the choice of the trial wave function. Sometimes as 
the number of parameters increases, the accuracy of WFO deteriorates, while the accuracy 
of energy is much more improved. Of course, in this case, the variational calculation becomes 
much more tedious and difficult. 


13 



The trial wave function with a single variational parameter is most convenient for use. 

4 

If the accuracy of 10% for WFO in the Cornell potential case is tolerable, ^ 

would be the best choice for the 15 state trial wave function. For the Martin and loga- 

3 

rithmic potentials, V’(^) = ^ is the most appropriate trial wave function for the 15 

state, and the corresponding WFOs have quite satisfactory accuracies. In particular, these 
forms can give very simple and reasonable relations between WFO and the reduced mass, 
which agree with those deduced from the general arguments for power-law potentials. How¬ 
ever, the requested accuracy of squared WFO in the system concerned is generally lower 
than 2 — 3%, and the trial wave functions with a single parameter cannot provide this accu¬ 
racy. For instance, one indicated in Ref.[^] that the estimated decay constant in Ref.[^], 
where a single-parameter harmonic oscillator trial wave function was employed in solving 
the Schrodinger equation in the Cornell potential case, may not be reliable, although the 
computed transitions and energy are quite reasonable. 


Here we would like to point out that to solve the Hamiltonian (Eq.(77) in Ref.[^]) 


= -1 a- 


m 


Cf Po o?g In rfi 
2 tt r 


(26) 


the exponential-type trial wave function with a single parameter (Eq.(79) in [1R|) 


h{r) = ^ e-T (27) 

02 

was employed in the variational framework. The paper reported that because the resultant 
energy agreed with the ’’accurate” value up to an order of ©(af), Eq.(26) could be an ap¬ 
propriate trial wave function. Then, the resultant WEO would be accurate enough to serve 
their main goal. However, our estimation shows that by using various parameters given in 
Eq.(26) (take = m in Eq.(27)), the relative deviations of squared WEO are about 0.40 and 
0.18 for cc and bh, respectively. These results would lead to improper theoretical predictions 
for the quantities which are closely related to WEO. 


As the conclusion we can draw, for the binding energy, most of trial wave functions, even 
with a single-parameter, can result a solution with a satisfactory accuracy. However, This 
is not always true for WEO. To obtain a reliable value of WEO, one has to adopt not only 
an appropriate trial wave function form for a specific potential, but also the proper number 
of variational parameters in the function. Our finding indicates that for a specific potential 
form and a not very higher accuracy of WEO, one can always find out a relatively simpler 
and more reliable trial wave function with an appropriate number of parameters. However, 
the form of the trial wave function strongly depends on the potential. In general, there is 
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no universal rule to determine the form and the number of the parameters of the trial wave 
function. 

On the other hand, our study shows that if one chooses Eq.(16) as the trial wave function, 
when the potential is flatter at the large r, the higher power terms, namely the higher 
configuration mixing, should be considered so that the higher accuracy of WFO can be 
reached. Moreover, the results by using Eqs.(17) and (18) indicate that if several components 
of the trial wave function can compromise the descriptions of the asymptotic behaviors of the 
potential at the short- and long-ranges, respectively, the trial wave function would have the 
simplest form and can provide higher accuracies for both binding energy and WFO. Namely, 
the trial wave function is more efficient. Usually, when a four-parameter trial wave function 
is chosen, in the commonly used potential cases, the accuracies of the energy and WFO can 
reach 10“^ ~ 10“® and 10“^ ~ 10“^, respectively. 

As mentioned in the introduction, the variational method should further be studied. 
Now, if the potential is not very singular, namely would not cause the divergence in solving 
Schrodinger equation, we find the way to construct an efficient trial wave function to get the 
more accurate binding energy and WFO. Then the next step is to study that if the potential 
is very singular, say more singular than 1/r^, how the bound state problem can be solved and 
the real physics can be obtained by using the variational method, and whether the obtained 
result in this method corresponds to reality. These discussion is shown in our next paper |^]. 
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